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The relativistic theories of light propagation are generalized by introducing two new 
parameters <; and rj in the second post-Newtonian (2PN) order, in addition to the pa¬ 
rameterized post-Newtonian parameters 7 and [3. This new 2PN parameterized (2PPN) 
formalism includes the non-stationary gravitational fields and the influences of all kinds 
of relativistic effects. The multipolar components of gravitating bodies are taken into 
account as well at the first post-Newtonian order. The equations of motion and their 
solutions of this 2PPN light propagation problem are obtained. Started from the defi¬ 
nition of a measurable quantity, a gauge-invariant angle between the directions of two 
incoming photons for a differential measurement in astrometric observation is discussed 
and its formula is derived. For a precision level of a few microacrsecond (/ias) for space 
astrometry missions in the near future, we further consider a model of angular measure¬ 
ment, LATOR-like missions. In this case, all terms with aimed at the accuracy of ~ 1/ias 
are estimated. 

Keywords: astrometry; reference systems; relativity. 
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1. Introduction 

With the development of observational techniques and the improvement of measure¬ 
ment methods, it is time for astrometry to unfold a new era indubitably. Presently, 
the satellite laser ranging like laser geodynamics satellite (LAGEOS) has achieved 
accuracies of about one millimeter,® the precision of lunar laser ranging (LLR) has 
approached one millimeter,® and the very long baseline interferometry (VLBI) has 
attained to the precision of 0.1 mas or even better.® Beyond the current stages, 
astrometric observation will be able to attain the accuracy of a few microarcsec- 
ond (/ras) or higher for some astrometric missions in the future, such the Gaia,® 
the Space Interferometry Mission (SIM),® the Square Kilometer Array (SKA),® the 
Laser Astrometric Test Of Relativity (LATOR)Pand the Beyond Einstein Advanced 
Goherent Optical Network (BEAGON).® This tendency requires a practical frame¬ 
work that can satisfy the accuracy of /ras. 
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In order to realize this purpose, there are three approaches: post-Newtonian 
(PN) method, post-Minkowskian (PM) method and Synge’s world function (SWF) 
method. Several authors have investigated the light propagation in the general rela¬ 
tivity (GR) and alternative theories of gravity by PN method, which assumes gravi¬ 
tational fields are weak and motions are slow. Richter & MatzneiP studied 2PN light 
deflection for one body in parameterized post-Newtonian (PPN) formalisirP^bJJ i^y 
introducing another parameter A at c~^ of gij , the space-space component of metric. 
HellingJi^ discussed the relativistic effects in astronomical timing measurements in 
the PPN formalism, which could be applied to VLBI, spacecraft ranging and pul¬ 
sar timing in the gravitational field of a motionless body. Brumber^i^ detailedly 
researched the 2PN light propagation of the Schwarzschild solution in three gauges 
(standard, harmonic and isotropic) through introducing two coordinate parame¬ 
ters. Brumber^i^ also studies the angle between two incoming rays in a motionless 
A^-body system in IPN approximation of GR. Klioner & KopeikirfSl developed a 
practical relativistic model for the 2PN light propagation in the harmonic gauge of 
GR, in which the IPN contributions from bodies in uniform motion and the 2PN 
contributions from the Sun were considered. Motivated by ESA’s Atomic Glock En¬ 
semble in Space (ACES)P^ on International Space Station (ISS), Blanchet et al.^^ 
discussed time delay and frequency shift for one body in the isotropic gauge at IPN 
of the framework of general relativity. Klionei^ presented a practical relativistic 
model for Gaia under IPN formalism with PPN parameters /3 and 7 for uniformly 
moving bodies. Xie & Huan^i^ deduced the 2PN approximation of Einstein-aether 
(7E) theory in the form of both superpotentials and an A^-point mass. Kopeikirfi^ 
derived an explicit Lorentz-invariant solution of the Einstein and null geodesic equa¬ 
tions for data processing of the time delay and ranging experiments for a moving 
body in PPN formalism by introducing a parameter <5 in spatial isotropic term of 
c~'^ for Qij. Minazzoli & ChauvinearP^ extended the IAU2000 resolutions to in¬ 
clude all the terms for the requirements from some space missions in GR. And 
then thejPi^ investigated scalar-tensor propagation of light in the inner Solar Sys¬ 
tem including relevant c~^ contributions for ranging and time transfer. Klioner 
& Zschock^^ formulated the light propagation in 2PN framework of a stationary 
gravitational field for the Schwarzschild metric (one body) in harmonic gauge by 
introducing one parameter ekz in spatial isotropic and anisotropic terms of for 
gij. Deng & Xi^^ investigated 2PN light propagation model at the c~^ level in the 
framework of the scalar-tensor (ST) theory for a gravitational A^-point mass system. 
And thejP^ found the parametrization with single parameter in spatial isotropic and 
anisotropic terms of for gij is not valid for the ST theory. 

Within PM method, it only assumes the condition of weak gravitational fields, 
but not necessarily slow motion, and it expanded quantities in powers of the gravi¬ 
tational constant G. Kopeikin & Schafei!^ derived the light propagation at the first 
PM (IPM) linearized Einstein equations in the gravitational field of an arbitrary- 
moving A^-body system. Kopeikin & MashhoorP^ researched the effects of bodies’ 
spins on the light propagation in the gravitational field of an arbitrary-moving 
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A^-body system in the harmonic gauge to 1PM order under GR. A model of the 
celestial sphere and of the observables for next-generation astrometric missions was 
constructed by numerical simulations with the 1PM method in GR for A^-body of 
the Solar System in two cases: static® and dynamical.® Kopeikin & Makaro\® 
took into account light deflection by a giant planet with multipoles in 1 /ras astro¬ 
nomical measurements using the 1PM method for two-body (the Sun and Jupiter 
or Saturn). 

In the method of SWF for the light propagation, the integration of geodesic 
equations can be avoided.®® Linet & Teyssandier® derived time transfer and 
frequency shift to the order c~^ in the field of an axisymmetric rotating body 
with this method. Le Poncin-Lafitte, Linet & Teyssandier® worked out a recursive 
procedure for expanding SWF method into a perturbative series of ascending powers 
of G. Le Poncin-Lahtte & Teyssandier® used SWF to investigate influence of mass 
multipolar moments on the deflection of a light ray by an isolated axisymmetric 
body. Teyssandier & Le Poncin-Lafitt^® derived the PM expansion of time transfer 
functions with SWF method. Recently, with the some method (SWF), Hees, Bertone 
& Le Poncin-Lafitt^® present a procedure to compute the relativistic coordinate 
time delay, Doppler and astrometric observables up to the second PM (2PM) order, 
and study light propagation in the field of a moving axisymmetric body for the Juno 
mission.® Bertone & Minazzoli et aP^ built the time of flight and tangent vectors 
in a closed form within the SWF formalism giving the case of a time-dependent 
metric and shown how to use this new approach to obtain a comparison of the time 
transfer with space astrometric modelings. 

In this work, we use the PN method. Table [T] gives a comparison among these 
previous works and our present work. Firstly, it is worth emphasizing that we employ 
the same integral technique and iterative method in Refs. [T31 [TU [T71 [HI Aimed at 
the level of /ras, our work has the following aspects: 

• N gravitating bodies in the Solar System, including the effects of their quadrupole 
moments and spins; 

• motions of the A^-body. Kopeikirf® pointed out and discussed that these mo¬ 
tions would contaminate the observed numerical value of PPN parameters. In 
our model, we separate them into two parts: the IPN contributions from the rec¬ 
tilinear and uniform A^-body (in a short time span) and the 2PN contributions 
from the motionless Sun, which is sufficient for near future missions; 

• two new 2PN parameters c and rj in addition to the two PPN parameters f3 and 

7 . and rj respectively parametrize the spatial isotropic and anisotropic terms in 
gij at In some cases, and rj are functions of /3 and 7 ,® but, in other 

cases, they can be independent of the PPN parameters.® 

In what follows, our conventions and notations generally follow those of Ref. |36l 
The plot of this paper is as follows. With the level of accuracy of ~ fias, we develop 
a 2PN parametrized (2PPN) framework for light propagation in section 2. Subse¬ 
quently, in section 3, a gauge-invariant angle between incoming light-rays is derived 
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and some special cases are shown. Then, we discuss a LATOR-like model of angle 
measurement: an observer measures the angle between two incoming rays emitted 
separately from two spacecrafts which are all at the Earth’s orbit circling the Sun. 
With the level of accuracy of ~ 1/ras, we estimate the magnitudes of contributions 
from various physical sources in this case. Finally, conclusion and perspectives are 
outlined in section 4. 

2. 2PPN framework for light propagation 
2.1. 2PPN metric for light 

In our investigation, we neglect the effects of our galaxy and external galaxies on 
the Solar System. This means the Solar System is an isolated system. Observations 
are modeled in the Solar System barycentric reference system (SSBRS), which can 
be mathematically described by a metric tensor. 

For practical reasons, we only care about the effects more than Igas in this 
work. It means the metric for light can be necessarily (but not over) simplified. The 
2PPN metric for SSBRS reads as 



where e = 1/c and 0{n) means of order e". mA and S\ are the mass and intrinsic 


angular momentum (spin) of the body A; subscript “ 0 ” denotes the evaluation 
of the Sun; is the coordinate velocity of the body A; r\ = — x\{t) and 

the trajectory of the body A is represented by x\(t) and ta = 
j<i/c> _ _ igikj^ jg symmetric trace-free (STF) quadrupole moment of 

body A. and rj are two new 2PPN parameters. 

At the IPN order, we keep all of the contributions from the mass monopoles of 
the Sun and the planets in the terms of goo and gij. Table [2] shows their effects 
of the major bodies in the Solar System on the light ray when the light ray grazes 
the limb of one body. The biggest one comes from the Sun, which is about 10“®. 
We also include the quadrupole moments at this order to take account their leading 
contributions. Table |3] shows the effect of non-spherical part (J 2 ) of a body on light 
when a light ray grazes its limb. Although the fourth (J 4 ) and the sixth (Jg) zonal 
spherical harmonics for the giant planets, such as Jupiter and Saturn, can cause 
deflections larger than 1 gas for light deflection when the light ray arrived their 
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limbs, ^ we leave them untouched in this work. In fact, their effects can be easily 
taken into account from the studies for Gaia mission.E^I^ 

In goi [see Eq. 0], we include all of the contributions from the mass monopoles 
and their velocities of the Sun and the planets in the terms of goi. Table [3] gives 
the estimation of got related to a body’s orbital motion v\. In order to investigate 
the influence of the gravitomagnetic fields on light propagation, not only a body’s 
orbital motion but also its spin are taken into account in the metric of goi in our 
model. 

For the 2PN terms in goo [see Eq. ([T])] and [see Eq. ([S])], we only consider the 
contributions of the monopole of the Sun and neglect the 2PN contributions due to 
the planets as well as the nonlinear combinations between the planets and the Sun, 
and the terms 0{e'^J2A) whose contributions are less than l/ras. In the e'^ terms of 
goo and gij, we also omit the terms explicitly depending on Vq, such as ^ 4 ™® Vq , 
because their contributions to the deflection are at the level of femto-arcsecond. 

It is very important to point out that the terms in goo and g^ and the 
terms in goi are all time-dependent, even including the terms associated with 
the Sun due to its barycentric motion Xq (t) and Vq (t) . KopeikiiP^ demonstrated 
that this motion of the Sun affects the measured values of the PPN parameters, and 
cannot be ignored. Moreover, the barycentric motion of the Sun was crucial in order 
to test general relativity in the Cassini experiment.!^ Thus, when we integrate the 
trajectory of the light-ray (see Sec. 2.3), we will keep the velocities of all the bodies 
including the Sun. 

At the 2PN order of g^ , we introduce two new parameters and g, which 
generalize the previous works. They have different values in different gravitational 
theories (see Table |5|). When <^ = g = 1 , the Eq. dS]) reduces to GR. When c = 
27 ^ — i -I- 2/3 — I and g = ^(1 -I- 7 ), the Eq. ([3]) returns to ST.!^ When = 1 -|- 
and g = 1 — ^Ci 4 , the Eq. (jS)) coincides with When <^ = g, the Eq. m goes 
back to the results of Ref. [22l When ^ = |(5 and rj = 0, the Eq. 0 coincides with 
the result of Ref. [191 


2.2. 2PPN equations of light 

Generally, for a photon propagating in a spacetime in which Einstein Equivalence 
Principle (EEP) is valid, the basic equations of lighfP^ are 


dx^^ dx'^ 

d^x^^ dx''dx^ 

-r''- 

dA2 "" dA dA 


= 0 , 

= 0 . 


(4) 

(5) 


We replace the affine parameter A with coordinate time t, and the equations 
become 


„ dx^^ da;'" 

0 = g/ii/ 


dt dt ’ 


( 6 ) 












May 26, 2015 0:7 WSPC/INSTRUCTION FILE IJMPDdeng20150523 


6 X.-M. Deng 




= -P*^ 


dt 


- r* 


dx'' dx'" 
dt dt 


( 7 ) 


Assuming x = cs/x = 0(1) and fx ■ fi = 1, we find the expression for s from Eq. 
©• Then, by substituting x • x = and the metric Eqs. m-m into Eq. © , we 
obtain the 2PPN equations of light propagation in SSBRS as follows 


X* — + J^Q + J-g d- 

where the IPN monopole component with the orbital motion is 


( 8 ) 


E® — — 
IPN — 




X ■ VA 


-2 1 - 


A ■ ^ 

X • vyi \ X • 






TCA ■ 


the influence of quadrupole moments of the bodies in the Solar System is 


^Q = 


3(1 + 7 ) G y<jk>^k 

A ^ 


-5 


r\r\ 


+ 2Sij-2{ 2—-5 


X ■ ta 


or I 


(9) 


CIO) 


the effect from their spins is 


G f (S^ 


X X 


= 2 ( 1 + 7 ) +^^[(Sa xr^)-x]r)t-^^[(SA xr^)-x](x-r^)x^ 


A ^A 


-(x-r^)(SA X r^)® L 


2 c^r 

and the 2PN monopole component of the Sun is 
G'^rn? 


( 11 ) 


T* — 9. 

-f 9PA7 — ^ 


"o 


'© 


- i(c + ry) + 27(1 + 7 ) + 




2(1 -/3) + ^ - 27^ 


(X^.. 


( 12 ) 


When we consider a case of stationary gravitational field for one body and <; = rj, 
our result Eq. (jS]) will be identical with Ref. [22] When we take GR into our model 
(7 = /3 = <j = 77 = 1), Eq. ([ 8 ]) will reduce to the result in Ref. [TH 


2.3. Trajectory of the light-ray in the 2PPN Framework 

The trajectory of the light-ray in the 2PPN framework can be obtained through 
integrating Eq. (| 8 |). However, it is difficult to derive it directly. So we adopt an 
iterative method used by Refs. [TSl [HI [22] 

But before that, we should find a way to describe the bodies’ motion analytically 
in the calculation. Chebyshev polynomial could perfectly deal with the motion of 
the bodies by interpolated ephemerides. This has been guaranteed at DE/LEj^ 
INPOF®! and EPIVP^ ephemeris , due to its convenience for numerical calculation. 
For analytical calculation, however, it is impractical. As an analytical ephemeris. 
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VSOpE^ used an iterative method by perturbations up to the the eighth order of 
masses for the planets in order to deal with the motion of the bodies. But this 
approximate method is still lengthy for our investigation. 

A Taylor expansion for the trajectory of the body A at some moment Ia is 
the simplest and practicable for our purpose. However, it leads to two sticky issues. 
One is the order of the Taylor expansion. The other is the determination of the 
moment tA- For the first issue, we suppose that the motion of the bodies, including 
the Sun, is rectilinear and uniform, that is, 

x^(<) = XA(tA) + v^(<a)(< - tA) + 0[{t - tA)^]. (13) 

Klioner & KopeikirfS! have pointed out that the residual terms of Eq. (IT3t are 
negligible for the accuracy of 1 /ras through a reasonable choice of the moment tA- 
This assumption, which means the path of a celestial body is approximated to a 
straight line, is enough for some measurements conducted in a short time span since 
the time of propagation of electromagnetic waves is very short with respect to the 
orbital period of a body in the Solar System. For instance, the time the light takes 
from the Sun to the Earth is about 8 minutes while the orbital period of the Earth 
is about 365 days. 

For the second issue, Klioner & KopeikiiP^l indicated that tA can be used to 
minimize the error in the solution of the light propagation equations caused by the 
higher-order terms neglected in Eq. (HI. Klioner & KopeikirfSlalso showed that the 
minimization procedure makes tA equal to the moment of the closest approach of 
the unperturbed light ray to the body deflecting the light ray. Kopeikin & Schafei!^ 
proved rigorously by solving Einstein and light ray propagation equations that at 
higher astrometric precision tA must be taken as the retarded instant of time cor¬ 
responding to the retarded (Lienard-Wiechert) solution of linearized gravitational 
field equations. Klioner & PeipP used the numerical simulations and showed that 
it was sufficient to use the well-known solution for the light propagation in the field 
of a motionless mass monopole for the accuracy of ~ 0.2 /ras and substituted in 
that solution the position of the body at the moment of closest approach. Klioner 
& PeijP found the post-Newtonian analytical solution for the body being at rest at 
its position at the moment of closest approach or at the retarded moment of time 
are virtually indistinguishable from each other for the Solar System applications 
and showed it attained an accuracy of 0.18 /ras if we took the form of uniformly 
moving bodies like Eq. (US. Therefore, we implement integration technique like 
Klioner & KopeikirP to derive the trajectory of light-ray. 

We assume the unperturbed light-ray as follows 

xat = xo -k c[t - to)n, (14) 

where fi is a unit vector representing the light direction at < = — oo, to is an instant 
on the light path and Xq is the position of photon at tp. The photon’s coordinates 
can be written as sum of perturbations with respect to xjy: 

x(t) = xat -k ^x = xat -k ^xipiv + (5 xq -k ^X 5 -k 6 x 2 pn- (15) 
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For J'lPTv, we use the following assumption for motion of the bodies 

YA{t) = - XA(t) = x(t) - XA(tA) - ■VAitA){t - Ia), (16) 

where Ia is the moment of the closest approach between the body A and the un¬ 
perturbed light ray. Due to the smallness of Tq, J-g and J^ 2 PN’ comparison with 
■^iPNi it i® sufficient to suppose that 

YAit) =XN{t) -XA{tA)- (17) 

After these, we can obtain the following results 

ix(t) = n -(5xiPAr(xAr -t (JxiPAf) -t -(5xQ(xAr) + -(5xs(xAr) -h i(5x2PAr(xAr0l8) 


x(t) = XAr(t) ■ 


(5xiPAr(xAr -I- (5xiPAr) - (5xiPAr(xo) 


+ 


(5xq(xjv) - (5xq(xo) 


(Ixs(xAr) - (5xs(xo) 




<5x2PAr(xjv) — 5x2PAr(xo) 


(19) 


It is worthy of note that the 2PN terms in our solution actually have two sources, 

direct and indirect. The direct part comes from the 2PN order itself. The indirect 

part comes from the IPN terms when the IPN solution is iterated into itself in order 

to attain a 2PN accuracy, namely, we substitute x^v + ^xipjv into the trajectory of 

the light-ray in the IPN approximation. During the procedure of iteration, we keep 

only the 2PN terms related to a motionless Sun. We can obtain that 

1 /fix (r^xk) 

-dxiPAr(x) =-(1-h 7) > ^- k' 

C ^ C^TA 


i(5xQ(x) = (1 + 7 )^ 


G 


A ^ 


kvA - k • TA 

rA(2r^ - n ■ r^) 
(r^ - n • vaY 


Pa + 


1-3 


(n • taY 




-3- ya -Sa 


ta 

i<5xs(x) = (1+7)^§| 

(Sa X n) 


Sa X dA <YA{2rA -n-rA) 
r\ r\(rA-'A-YAY 


(Sa X n) - ^ 
a A 


( 20 ) 


( 21 ) 


Sa • (n X dA) 


taYa - n- rA) 


( 22 ) 


1 .. , , 1 GWq 

-6x2pnA) = --T] 4^4 (n • r0)r0 


4 4 


G^ml 

- -A -<4(7) 


(1 + 7)" 




r©(r0 - fi • rg) Vrg r© - n • r© 


+ [/3 - 7 ( 2 ^ + »7) - 2(1 + 7)]-)2- 


-^0 


n • r, 


© 


'© 


(Atq 


-n 


1 / TT n • r© 

— — -I- arctan —-— 

d© \2 d© 


1 


27(1 + 7) +/3 - 7*^) — - (1 + 7)^- 

2 yr© r©-n-r©j 


(23) 


and 

Sxipjv(x) = -(1 + 7) 


E Gttia r ^4 

1 Ir 


X (rA X k) 


\ krA — k • 


+ k In kvA + k • 


(24) 
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( 5 xq(x) = (1 + 7 )^ 


G 




+ n ■ 

rA-n-VA' 


Pa 


d\{h-rA) 


d\ 

<\A + -T'^A - 


n ■ VA 

ta 




^xs(x) = (1 + 7 )^ 


G (n-VA 


+ - 


A 

n • ta 


(^ta 


d\ 


TA-n-VA 


(Sa X n) 


(5x2PAr(x) = -rj 


1 G^m: 


2^2 

0 


A'' 


-(1 + 7)^ 


r© 


G^mg 


(Sa X dA) - 


1 


dA(r-A + n ■ ta) 
d\{rA - n • ta) 


[Sa • (n X dA)] 


(26) 


/3 - -(2^ + 77 ) - 2(1 + 7 ) — arctan 


dpi 


n ■ rg 

dg 


rg -n-rg 


G^rri^ 


d© ( 1 + 7 )^ 


1 


(r© -n-rg)2 


1 


+ [/^ “ ^( 2 *? + 7 ) - 2(1 + 7 )]—^ 


arctan 


dp 


(27) 


where dA = n x {ta x n) is an impact parameter for the body A which is the closest 
approach of the unperturbed light ray, dA = |dA|, ta = x^v — XA(t), va = Ita], 
k = n — VA(tA)/c, k = |k| and 


p), = 2 jf 


r<3k> 


Wd\ 


_ O T 

c\a — 


d\ + r. 


■<La 

ta 


P, = 


nJf^k -^A 

d\ 




d\ 


A 
rA '' 


s\ = 2J, 


<^k>f^k _ ^j<jk> 


Wdh 


^d\ - 


.2 ^A 

“a 




d^A 


— 2 


d^j^d'X 

d\ 


Ia 

rA '' 


(28) 

(29) 

(30) 

(31) 


If we only consider a static case for one body and assume 77 = c, Eqs. (ESI) 
and EH) will return to the results of Ref. |22l When we neglect and Tg, 

Eqs. ESjl and E^ go back to the results of Ref. [ITl When we consider GR 
(7 = /3 = <7 = 77 = 1), Eqs. (fT^ and (IT^ will reduce to the results of Ref. [T4| 


3. Angular measurement 

3.1. Construction of a gauge-invariant angle in the 2PPN 
framework 

In practical astronomic measurements, a differential measurement is a more pow¬ 
erful method. This concept is employed by LATOR missiorP through a skinny tri¬ 
angle formed by two spacecrafts and the ISS. In what follows, we mainly focus 
on discussing LATOR-like missions. Firstly, we construct a gauge-invariant angle 6 
between the directions of the two incoming photons based on Ref. [131 It reads 


cos 9 = 


(32) 
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where the spatial projection operator is 

ha(3 — 9a(3 

which projects the two incoming photons onto the hypersurface orthogonal to the 
observer’s four-velocity u°‘ = dx^/cdr and iFf = dxf{t)/dt, and = dx 2 (t)/dt 
are the tangent vectors of the paths a;“(t) and x^it) of the two incoming photons. 
Then, we obtain 

(1) /(2) (2) (2) \ /(3) (3) (3) \ 

COsO = (fii • 62 ) + f obs + if obs + f IPN + f Q j + I f obs + f obsxlPN + / S 1 

/(4) (4) (4) \ 

d” ( / obs f obsxlPN f 2PN 1 ? 

where the first term is the angle between the unperturbed paths of the photons 

(n) 

from two signals, / denotes terms of order e", the subscript “o 6 s” denotes terms 
related to the observer’s velocity, subscript “IPN” denotes the contribution from 
the monopoles and orbital motions of gravitating bodies, subscript “Q” denotes 
the terms from their quadrupole moments, subscript “obs x IPN” denotes the 
coupling terms of the bodies’ mass and observer’s velocity, subscript “S” denotes 
the terms of bodies’ spin and subscript “ 0 ” denotes the Sun monopole contribution. 
The expressions for these functions are presented in [Appendix A[ Obviously, the 
position of the photon at the moment t of observation coincides with the position 
of the spacecraft (observer) so that 

Xobs = xoi -I- c{t - toi)ni + <5xi = Xo 2 + c{t - ^ 02)62 + ( 5 x 2 , (35) 

where (^ 017 X 91 ) denotes the moment and position of the light signal 1 of emission 
and (to 2 ,xo 2 ) for the light signal 2 respectively. Eq. (1^ also give a constraint at 
the moment of measurement that 


^lA — ^2A — ^obsA — 

Eq. (l34l) can be expanded with respect to the parameter e as 

(1) /(2) (2) (2) \ /(3) (3) (3) \ 

6{t) = obs + ( "d obs + "d IPN + "d Q j + (^?obs+ "d OM + "0 S j 

(4) (4) \ 

obs + "0 2PN 7 


(36) 


(37) 


where G ( 07 '^) is the angle between the unperturbed light paths from two given 
sources 


■00 = arccos(rii • 112 ), (38) 

(n) 

and 10 obs is the deflection angle due to the observer’s motion in terms of order 
( 2 ) 

£"■7 is the IPN deflection angle due to the spherically symmetric field of 
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(2) (3) 

gravitating bodies, "d q is due to quadrupole moment, om is due to the orbital 
(3) (4) 

motions of the bodies, s is due to the spin of the bodies and e"* i9 2 PN is the 2PN 

deflection angle due to the spherically symmetric field of the Sun (see [Appendix B| 
for detailed expressions). 

3.2. Special cases 

As one kind of the verification, our results could reduce to some special cases and 
some known results, such as aberration of special relativity, IPN light deflection, 

2PN light deflection and Light deflection due to spins. 


3.2.1. Aberration 

In Eq. ((571) . we pick up the terms due to the velocity of the observer 

( 1 ) ( 2 ) 

9 = Aq + 9 obs + 9 obst (39) 

and assume that one signal propagates along the direction of motion of the observer, 
namely, Vobs = ffii. Then, we obtain 

9 = Aq + ev sin'do + sini?o cosAq. (40) 

This is just the aberration given by special relativity after ignoring 


3.2.2. IPN light deflection 

The path of light is bent in a gravitational field, which represents by 

( 2 ) 

9=9 IPN- (41) 

We assume that only the Sun is considered and one of the two light rays is just along 
the line connecting the Sun and the observer and without any bending, namely 
source 2. Then, we obtain 

(/1 + 7 \ 4 Gm 0 /1 + cos 9o 

where we use n 2 • di^/diA = sin do, rii • vaIta — cosdo and = t'a~ (^i ' i’a)^- 
For GR (7 = 1 ), When cos-do = 1, our results reduces to the famous infinity-infinity 
light deflection formula. 

Brumber^i^ discussed the angle of the two incoming photons in IPN general 
relativity for a A^-body system. Now, we discuss our result about IPN gravitational 
deflection, namely, 

( 2 ) 

cosd = (ill • n 2 ) -b / IPN 
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= (fii • 112 ) - e^(l + 7 ) 
= (fii •fi2) + e^(l + 7) 
When 7 = 1 , this is just the 


E 

A 

E 


GruA ( ni • [n 2 x (r^ x n 2 )] 62 • [fii x (r^ x fii)] 


ta 

GniA 

ta 


\ 


ta - n2 • VA 
ta X fii ta X fi 2 


rA - ni ■ ya ta - n2 • ya 


rA - ni • Ya 

(fii X 62 ). (43) 


result of Ref. [131 


3.2.3. 2PN light deflection 

By extending the previous calculation to 2PN order, we have arctan(ni -YQ/dio) ~ 
7r/2, diQ/rQ Ri 0 and 

( 4 ) r 1 I (72^2 

^? 2 PJV « I - 2(1 + 7)2 + [ 2(1 + 7 ) -/3 + -( 2 c + 7 )] 7 rj ^ 4 ^ 2 ® - (44) 

If the above result returns to GR = (3 = <; = rj = 1), it coincides with the result 
of GR within harmonic gauge . If (2? + 77 )= 3A, it will coincide with the result of 
Ref. m In the case of c = 77, it will reproduce the result of Ref. |22l 


3.2.4. Light deflection due to spins 

The spin of a body can also produce a light deflection. It belongs to one part of 
gravitomagnetic field. It reads as 

(3) 

9=ds. (45) 


Repeating previous calculations, we obtain 


e Ri 2(1 + 7 ) 


„3a2 

C 


+ 4(1 + 7 ) 


[Sa ■ (ni X diA)] 
C^d^A 


(46) 


This is just the result of Ref. [25| with 7=1 with turning their vector angle for 
deflection by reason of spin to a scalar quantity. 


3.3. A MODEL FOR LATOR-LIKE MISSIONS 

In this section, we apply our model to LATOR-like missions by qualitative esti¬ 
mate.We assume there are three spacecrafts in an orbit circling the Sun at the same 
distance as the Earth. Two of them carry a light signal emitter on board and the 
third one carries a light signal receiver. When the two spacecrafts and the third one 
are on the opposite side of the Sun, the experiment will be conducted. When a light 
signal 1 emitted by one of two spacecrafts passes by the limb of the Sun, another 
signal 2 is emitted by the other spacecraft which has quite a distance from the Sun. 
The receiver located on the third spacecraft which is on the opposite side of those 
two spacecrafts measures the angle between the two incoming photons. In the case 
of LATORj^the observer will be set on the ISS. 
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We assume that the velocity of the observer in SSBRS is near the orbital velocity 
of the Earth. In this case, we need to consider three bodies’ gravitational fields: the 
Sun, Mercury and Venus. We can neglect the effect of gravitational field of the 
Earth through making the position of the observer away from the Earth. However, 
for the case of ISS, the gravitational field of the Earth affects the motion of ISS and 
this effect can be represented by equations of motion for ISS. The influence of other 
gravitating bodies in the Solar System on the light propagation in this case can be 
neglected. For LATOR, signals are emitted inside the Solar System, but we could 
extended the trajectories to t = —oo to make make ni and 112 meaningful. With 
the angle between the initial emitting directions of two light signals 1 and 2 


= fii • n2 R::! 1° 


(47) 


it means that the distance between these two spacecrafts is about 5.22 x 10® m. We 
can simplify Eq. dSzl) for this practical case as follows 




= ■do + e(n 2 • Vof,^) tan — - e tan — + e 2(1 + 7 ) 


dr 


Gulp 


■0 


diQ 


T erml 


Term2 


+ e^2{l + 7 )^ + £ 22(1 + + £ 22(1 + 7 ) 


-^20 


^IMer 


Term3 

Gray 

div 


TermA 


Termb 


+ £ 22(1 + 7 )^^ - e" 2 (l + 7 ) tan 
d2V 


Term6 

doN GniQ Vo5s ■ di0 


Ji 0 


d\p 


10 


Term! 


Terms 

2 ^ 


+ £^2(1 + 7)—; —^(n 2 • Vobs) — £“^2(1 + 7 ) 
« 1 © 


'•0 


di 0 


Term9 TermlO 

- - 7(2‘r + 1?) - 2(1 + 7 )^]^^ 

4 «!© 


0 


—h arctan 


ni ■ Tq 

diQ 


1 


Termll 

m, 


-£"[/?-02,+ ,y)-2(1+ 7 )^]^^ 

4 “20 


—h arctan 


n 2 ■ r© 

d2Q 


Terml2 


+0{< l/itts). 


(48) 


with the cut-off precision of ^ 1 /las. Here we consider the largest influence of 
Mercury (subscript “Mer”) and Venus (subscript “V”) on the measurement when 
the light signals might pass by the limbs of them. We estimate the order of these 
terms on Eq. (l 4 ^ which are listed on Table [6l 

For LATOR mission, it will attain the precision of 0.01 /ras. It needs to con¬ 
sider a complete form at c~'^ of gij and goQ. Various terms at of gij and goo 


respectively are 


GrriA „.2 

— 3 — u 

C^VA 


A’ 


E. 


EaE 


B^A 


GmAmB 

C^rATAB ■ 




B^A 


GrriAmB 

C^rATB ' 
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E GrriA 

A c*r\ A' A' 




_ GruAmB _ „3 

B^A c‘'‘r\^(rA+rB+rAB) ^B 


GniAms 


-_l—ay These terms have be shown 

' TAB TA ' ' tab rB ' 


EaEb^i^A c*{rA+rB+rAB 
by Ref. 1231 These terms can be considered as different combinations between terms 

in Tableland in Table 01 If the accuracy of the model is required to reach 0.01 
fj,as, we need to deal with two problems: ( 1 ) integrating the two-body terms in the 
2PN order and (2) introducing new parameters of the coupling terms in the 2PN 
order. Besides, LATOR also measures laser ranging for three sides of the triangle, 
it needs to derive the light time solution. And the spacecrafts in LATOR mission 
are constantly moving. The triangle formed by the spacecrafts are also changing. 

~ 1 ° is just a special case which presents the maximal effect of this measurement. 
We will study these issues numerically in our next move. 


4. CONCLUSIONS AND PROSPECTS 

In this paper, we present light propagation and a gauge-invariant angular measure¬ 
ment in the 2PPN framework by introducing two new parameters c and rj besides 
the two PPN parameters 7 and (3. In the framework, we consider all kinds of rela¬ 
tivistic effects on light propagation in SSBRS, which are monopole and quadrupole 
moments of the bodies in the Solar System, their motions and their gravitomagnetic 
fields. With the derivation of a gauge-invariant angle between the directions of two 
incoming photons, we further discuss a practical astronomic observation, namely, 
an observer on a spacecraft measures the angle between the two incoming photons 
emitted separately two spacecrafts which are all at an orbit circling the Sun at the 
same distance as the Earth. Given attaining a level of ~ 1 /ras for space astrometry 
missions in the near future, like LATOR mission, the terms at this level are listed. 

In this work. We approximate planetary motions as linear motion. This approxi¬ 
mation is not valid for the time scale comparable with the orbital periods of planets. 
Our next move is to study this problem numerically by integrating the equations of 
motion of A-body and light simultaneously. How to parameterize the 2PN coupling 
terms in these experiments is another issue that will be investigated. 
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Appendix A. Functions in cos 6 

The expressions of the functions in Eq. (1341) are 

(fii • n2) - 1 


( 1 ) 


f obs = 


(A.l) 
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( 2 ) 


• Vobs)^ + 9 ^ ■ Vo&s)(nj • -Vobs) - 


k=l k,j=l k^j 

( 2 ) ^ ^ 1 


vlbs I I (ni • n2) (Ai: 


/ IPN — EE hj ■ -(5xiPAr(xAr + 5xipjv)|^ 


h,j=l h^j 
2 


— ^(fii • n2)nj • -i5xiPAr(xAr + 5xiPAr)|j, 


(A.3) 


( 2 ) ^ ^ 1 1 


■^<3= E E"j ■ “E("1 ■ "2)nj •-'5 xq(xjv)|j, (A.4) 


(3) 

fobs=^ 


h,j=l h^j 

2 2 


i=i 


E(“fc " - E("fc ■ '^obs)l’obs A E E("fc ■ ^obs) (^j * ^obs] 


k=l 


k=l 


k,j=l k^j 


(fil • n2) - 1 


(A.5) 


(3) ,A 1 . 

/ obsxlPN = ^N^obs • -^XiP7v(xjv)|fc 
fc=l 


(fil • fi2) - 1 


+ e ^ (fi; • Vobs)nj ■ -(5xiPAr(xjv)U 


j,h,l=l 


2 j 

-3e^(fii • fi2)(fij • Vof,s)fij • -(5xiPAr(xAr)|j 


J =1 

2 


- E E("3 • fi2)(fii • Vobs)fi/j • -(5xiPAr(xjv)U, 

j,h=lj^h ^ 


(A.6) 


(3) ^ ^ 

fs= E E"j ■ “E(“1 ■ "2)nj •-<5xs(xiv)|i, (A.7) 


(4) 

f obs=^ 


j,h=l j^h 

2 2 


i=i 


^(fifc • Vobs)-* - ^(fifc • Vobs)^vL - 9 E E^"'' ■ ^obs){nj ■ Vobs) 


-^obs 


2 






+ E E("^ * (lij • Vq^s) H- 9 E E("fc • Vobs)^(fij • Vobs)" 


k,j=l k^j 


k,j=i k^j 


(fil • fi2) (A18) 


/ ObsxlPN = -(1 + 7)e‘‘ E ■ Vo&^)^(ni • n 2 ) + (1 + l)t^ ^ ^ ^^(fifc • Vobs)(ni ■ Vobs) 




't'obsQ 




. ^obsQ 


2 1 2 ^ 

-e^ ^ • -i5xiPAr(xAr)|b + 3e^ E(”l ■ ’^o6s)(fi2 ' ^obs)^j ' -(5xiPAr( xat)| j 


3,h=l 


i=i 
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2 1 2 ^ 

+3e^^(nj ■Vobs)'^iij ■ -SxiPN{xN)\j + X! ■ ^obs)'^nj ■ -6xiPN{xN)\h 

i=l ^ j,h,l=lj^h 

2 I 

+e^ X!(ni • Xobsfi^h ■ -6xiPN{xN)\h 

j,h=lj^h 

2 j 

+e^ y^(ni • Vo&s)(n 2 • Xobs)nj ■ -SxiPNixN)\h 

j,h=lj^h ^ 

2 1 2 

-2e2 y^(nj • Vobs)'Vobs ■ -SxipN{xN)\j - X! ' '^obs)Xobs ■ -S±iPNi^N)\h 

i=l j,h=lj^h ^ 

2 1 2 

+26^ E i;obg(ni • n2)nj • -5xipjv(xAr)|j + 2e^ E("i • n2)(nj • Vobs)vobs ■ -6xiPNixN)\j 
1=1 ^ 1=1 ^ 

2 I 

+e^ E E^^i ■ ”2)(nj • Vobs)xobs ■ -SxiPN{xN)\h 


3,h=lj^h 

2 1 

E E(-i • n2)(nj • Vobsf^h ■ -5xiPN{y^N)\h 

j,h=lj^h ^ 

2 j 

-3e^y^(ni • Vobs)(n2 ■ Vobs)(ni ■ n2)nj ■ -Sxipjv(xjv)lj 
1=1 ^ 

2 ^ 

-6e^E(*i ■ ft2)(ni • Vo6s)^fij • -(5xiPAr(xAr)|j, 

1=1 ^ 


(A.9) 


f2PN = 2 ^e 2 ^ 2 ^ 4 (”fc ■ i’ohs©)(nj ■ r^bsQ) - 2 ^ 4 K ’ >'ob«o) ("i ' 112) 

k,j=l k^j ^obsQ ^o6s© 


2 1 

- E E-i- -(5xiPAr(xAr)|jni • -5xiPAr(xAr)|^ 


j,l,n—l l^n 
2 


1 . ^ 1 1 
—-E(*ii • n2)-(5xiPAr(x]v)|j • -(IxiPAr(xAr)lj 
^ 1=1 
2 

3 , ^ 1 111 

+ 9E("1 • n2)nj • -(5xiPAr(x7v)|jnj • -(IxiPAr(xAr)lj + -(5xiPAr(xAr)|i • -(5xiPAr(xAr)|2 


1=1 


1 1 ^ ^ ^ 

+ (ni • n2)n2 • -(IxiPAr(xAr)|2ni • -'JxiPAr(xjv)|i + EEni- -(5x2PAr(xAr)U 


3,h=lj^h 


-y^(ni • n 2 )nj • -5 x2PJv(xjv)|4, 


(A.IO) 
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where C(x)|i and C(x )|2 denote that this function C(x) takes respectively the value 
at the signal 1 and the signal 2 . 


Appendix B. Angle expressions 

We list the expressions of the deflection angles in Eq. (EZl). 

^ , f- ,, 1 -cos-do 

d obs = e[(ni • Vobs) + (n2 • Vob^jJ 


d oh., = 


sin do 


sin do 

(ni • + (n 2 • ) ( 1 + ^ cosdc 


(B.l) 


(1 + cosdo)^ 

+ (ni • Vobs)(n 2 ■ Vobs) - Vobs(^ + cosdo) 


1 


(B.2) 


= (1 + y [ fti - [1^2 X (r^ X liOl ^ n.-|,iix (rpxni( 

^ TA Sin do t - n 2 • “ " " 


( 2 ) _ 

Q = -(1 + 7 )e" ^ ,3 \ K - (ni ' 62)^^ 


- ni • VA 

2rA -n2-rA 


o J<jk> 

■^:^(fii-d2x4) 

2 ri 


2n^d2x4 + (*2 ■ rA)fiin2 - {^2 ■ ta) 


{ta - n 2 ■ r^) 

diid^ 


2 “ 2 x 4 “2 


'‘2A“2x4 


(fi2 • di^) 


^2A J 

2h{d^A + (ni • rA)n{h^ - (ni ■ 


^Ia J 


- (fii • n 2 )ni] 


2rA-n,-VA k pfc 
/ - 

(rA-ni’VAY 


(B.4) 


(3) 

d obs — € 


3(1 + 2 cos^ do){l — cos do)^ 


+e- 


6 sin® do 

;COsdo(l — cosdo)^ 


n 3 


(fii • ^^obs) + (n2 • Vobs) 


sin® do 


[(fii • ^obs) + (fi2 -^obs)] 


vlbs - (fii ■ ^obsf - (fi 2 • Vobs)^ 


-(fii • Vobs)(n 2 • Vobs) +e® 


1 — cos do 


(fii • Vobsf + (fi2 • Vobsf 


sin do 

-(fii • Vobs)vli,s - (fi2 • ^obs)vlhs + (fii ■ Vobs)^(fi2 ' '•Tobs) + (fi2 ’ Vo{,s)^(fii • Vobs) 


+ (1 + 7)e' 


3 1 — COS do GniA 

sin do ^ ta 


(fii • Vobs) + (fi2 • Vobs) - 


'fobs ■ diA 'fobs ■ d2A 


+ (1 + 7)e^3y^ X! —^[(fii ■ Wbs) + (fi2 • Vobs)] 


sin do ^ ta 


TA-ni- Ta ta - n 2 • 

fii • d2A , fi2 • diA 


(B.5) 


(3) 


'doM = (1 +7)e^X! “ 


GlTlA 


sin do 


rA-n2-rA r^-nirA, 

(fii ■ d2A)(fi2 ■ va) _ (fi2 ■ diA)(fii ■ va) (fii ■ d2A)(fi2 ' va) 

VA^rA - n2 ■ Ta) rA(rA - fii • r a) (rA - fi2 ■ ta)® 
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^ (n 2 ■ diA)(ni ■ wa) _ (ni ■ d 2 A)(rA ■ v^) _ (n 2 ■ diA)(rA ■ va) 

(r^-ni-r^)2 r^(r^ - fia ■ r^)2 r^(rA - ni • r^)2 

(ni ■r^)(n2 ■ va) _ (62 ■ r^)(ni ■ v^) (n2 ■ r^)(ni ■ v^) (fli ■r^)(n2 ■ va) 
rAirA-m-YA) TAirA-'O-i-YA) rA{rA-'n2-YA) TAirA-ii-l-VA) 

02 -v^ - ,ni-VA . sn 2 -VAl 

+ (ni-n 2 )-|-(ni-n 2 )- , (B. 6 ) 


ta 


ta 


ya 


ta 


(3) 


'd s = -(1 +7)e^X! 


sin i9o 


f Sa • (fii X 02)] 

\ 1 (fi2-r^)' 

’(fi2 • Ya) (fi2 ■ Ya)'^ , 

9 ~r J- 

1 d2A 

YA 

YA Y\ 


+ 


fii • (Sa X d2A) Sa ■ (n2 x ni)] 


^Ia 


1 + 


(ni • r^) 


ta 


(ni • r^) (fii • 


+ 1 


ta 


02 ■ di^ 

dtA 

ni ■ d2A 

d2A 


[S^ ■ (oi X di^)] 


[S^ ■ (02 X d 2 A)] 


2 - 

2 - 


Oi • ya 


ya 

n 2 • Ya 


ya 


l + 2 M£d+("- , 7 )' 

ya 

ya Yi 


02 • (Sa X diA) .. oi ■ (Sa X diA) .. „ 02 ■ (Sa X d2A) 1 

+- 3 - “1 • ^2) - 3 -ni • 02 -3- 

t^A n J 


(B.7) 


4cos-do(3 + 2cos2z?o)(l-cos-do)‘‘p,. . .,4 

'dobs = -e -^ . 7 ; -^[(ni -Vobs) + (n2 • 


+e' 


8 sin’^ do 
4 (1 + 2 cos^ ^o)(l — cosdo)^ 


2 sin^ do 

+ (oi • Vohs)(o2 ■ Vobs) - Uohs] - e 


[(fii • Vobs) + (02 • Vobs)]^[(oi • Vobs)'^ + (62 • Vobs)'^ 


2 4COSl?o(l — COSI^o)^ 


sin^ do 


[(oi • Vobs) + (62 • Vobs)] 


x[(oi • Vobs)^ + (62 ■ Vobs)^ - (fii • Vobs)v'l - (02 • Vobs)vlbs + (^1 ' Vo6s)^(fi2 ' Vobs) 

/- n2/- m 4COSz9o( 1 - COSt?o)^ r/- ^2 ^2 

+ (n2 • Vobs) (ni ■ Vobs)] - e - ^ . 3 . -^[(ni • Vobs) + {^2 ■ Vots) 


+ (fii • Vobs)(n2 ■ Vobs) - Vobs] + e 


2 sin^ do 

2 i2 , 4I - cos-do 


[(fii • Vobs)"* + (n 2 • Vobs)"^ 


sin-do 

-(fii • Vobs)^Wobs - (^2 • ^obs)‘^vlbs - (fii • Vobs)(fi2 ' ^ obs)vlbs + («! ' ^ obs)^ {^2 ' Vobs) 
+ (62 • Vobs)^(ni • Vobs) + (ni • Vobs)^(n2 • Vobs)^] 




sin-^ do Yq 

'^obs * cil 0 Vobs ■ CI 20 


(fii • Vobs) + (62 • Vobs) 


4 (1 + 2cos^ do — 3 cos do) (1 — cos do) Guiq 


yq-ui-Yq r 0 -O 2 -r 0 j 

x[(fii • Vobs) + (*2 ■ Vobs)]^ 


+ (1 + 7 )e' 


fii • d20 02 • di0 


2 sin® do yq 

4 (cos do — 1 ) cos do GniQ 


Lr0-O2-r0 rQ-Oi-r0j 
X [(fii ■ Vobs)^ + (n2 ■ Vobs)^ + (fii • Vobs)(fi 2 • Vobs) - Vobs] 


+ (1+7K 

fii • d 20 _ 

JQ -h2-YQ^ rQ-h^■YQ 


sin® do Yq 

fi 2 • di 0 
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GmQ 


+ (1 + 7 )e'‘ cot '& 0 - 


re 


- 2(ni • Vobs) - 2 (n 2 • Vobs) - 2(ni • Vobs)in 2 ■ Vobs) 


^ (n2 ■ Vobg)(Vobg ■ dig) (til ■ Vobs)(Vofcs ' d20) 2 *-”2 ' '^of>s){'^obs ' d20) 


+2 


r© - Ill • r© 

(til • Vobs){^obs ■ di©) 


r© - 112 • r© 


r© - 112 • r© 


r© - ni ■ r© 


1 C^mr> 

+ (1 +7)e'‘^— -^[(fii • Vo6s)(n2 • ^obs) - v^bs + (ni • ^obs)'^ 


+ (n2 • Vobs)^] ( 2 + 


Ill • d 


2 © 


sin-do f© 

n 2 • di© 


r© - 112 • r© r© - ni • r© 


^ (fli ■ Vob^)(Vob^ ■ di©) ^ (n2 ■ Vobs){^obs ■ d2©) 


sin §0 r© 

^ (fii ■ Vobs)i^obs ■ d2©) (n 2 ■ Vobs)(Vohs • di©) 


r© - ni • r© 


r© - n 2 • r© 


r© - 112 • r© 


r© - ni • r© 


(B. 8 ) 


1? 2PN = —(1 + ^ 

sin 1/0 


n 2 • di© 

“i© 


1 + 


■ 2 © 


-[/3--(2^ + 7)+7"-l]e" 


fli • !•© \ fli ■ d 

^0 / *^ 2 © 

(ni • !■©) , (n 2 -r©)' 


1 + 


ii 2 ■ r© 


'© 


-[/5-4(22 + 7)-2(l + 7)]e 


m\ 


■0 


c?l© 

ni • da© 


C?2© 


+ 


n2 • d 


1 © 


'^ 1 © 


arctan 


sin i9o I di© 
ni_r© 

^ 1 © 


■ 2 © 


— + arctan 


^2 -i-Q 

d 2 © 


(B.9) 
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Table 1. Differences among the metric between some previous works and our work for light prop¬ 
agation. These include the order of expansion, whether or not to be parametrized, whether or not 
to consider spin, quadrupole moment, A^-body and bodies’ motion. 


Ref. 

Parametrized 

Parameter 

in 

Parameter 
— 4 

me 

Af-body 

Bodies’ 

motions 

Quadrupole 

moment 

Spin 

m 

Yes 

7 

/3,A 

No 

No 

Yes 

Yes 

EU 

Yes 

7 

0 

No 

No 

No 

No 

[m 

No 


— 

No 

No 

No 

No 

[m 

No 

— 

— 

Yes 

No 

No 

No 

RH 

No 

— 

— 

Yes 

Yes 

Yes 

Yes 

[Tel 

No 

— 

— 

No 

No 

No 

No 

nn 

Yes 

7 

0 

Yes 

Yes 

Yes 

No 

[Tsl 

No 

— 

Ci 4 

Yes 

No 

No 

No 

[THl 

Yes 

7 

0,6 

No 

Yes 

No 

No 

[ 20 ] 

No 



No 

No 

No 

No 

[U 

No 

7 

1,0 

No 

No 

No 

No 

EU 

Yes 

7 

/^) ^kz 

No 

No 

No 

No 

EU 

No 

7 

1,0 

Yes 

No 

No 

No 

EH 

No 

— 

— 

Yes 

Yes 

No 

No 

EH 

No 

— 

— 

Yes 

Yes 

No 

Yes 

EH 

No 

— 

— 

No 

Yes 

Yes 

No 

EH 

Yes 

7 

— 

No 

No 

Yes 

Yes 

EH 

Yes 

7 

0,6 

No 

No 

No 

No 

ED 

Yes 

7 

— 

No 

No 

Yes 

No 

[3^ 

Yes 

7 

0,6 

No 

No 

No 

No 

EH 

Yes 

7 

/5) ^hbl 

No 

No 

No 

No 

EH 

No 

— 

— 

No 

Yes 

Yes 

No 

EH 

Yes 

7 

— 

Yes 

Yes 

No 

No 

Our work 

Yes 

7 

0,';,v 

Yes 

Yes 

Yes 

Yes 


Table 2. The order of between bodies in the Solar System. This table shows that the bodies at columns affect the bodies at rows when 

c^vab 

a light ray arrives their limbs. 


GrriA 

Sun 

Mercury 

Venus 

Earth 

Mars 

Jupiter 

Saturn 

Uranus 

Neptune 

Sun 

2.1 X 10“® 

2.6 X 10-® 

1.4 X 10-® 

9.9 X 10-® 

6.5 X 10-® 

1.9 X 10-® 

1.0 X 10-® 

5.1 X 10-1® 

3.3 X 10-1® 

Mercury 

4.2 X IQ-i® 

1.0 X 10-1® 

4.9 X 10-1® 

2.7 X 10-1® 

1.4 X 10-1® 

3.4 X 10-1° 

1.8 X 10-1° 

8.7 X 10-11 

5.5 X 10-11 

Venus 

3.3 X 10-1'^ 

7.2 X 10-11 

6.0 X 10-1® 

8.7 X 10-11 

3.0 X 10-11 

5.4 X 10-1® 

2.7 X 10-1® 

1.3 X 10-1® 

8.2 X 10-1° 

Earth 

3.0 X 

4.8 X 10-11 

1.1 X 10-1® 

7.0 X 10-1° 

5.7 X 10-11 

7.1 X 10-1® 

3.5 X 10-1® 

1.6 X 10-1® 

1.0 X 10-1® 

Mars 

2.1 X 10-1® 

2.8 X 10-1® 

4.0 X 10-1® 

6.1 X 10-1® 

1.4 X 10-1° 

8.7 X 10-1® 

4.0 X 10-1® 

1.8 X 10-1° 

1.1 X 10-1° 

Jupiter 

1.8 X 10-12 

2.0 X 10-12 

2.1 X 10-12 

2.2 X 10-12 

2.6 X 10-12 

2.0 X 10-® 

2.2 X 10-12 

6.7 X 10-1® 

3.8 X 10-1® 

Saturn 

3.0 X 10-1® 

3.1 X 10-1® 

3.2 X 10-1® 

3.3 X 10-1® 

3.5 X 10-1® 

6.5 X 10-1® 

7.0 X 10-® 

2.9 X 10-1® 

1.4 X 10-1® 

Uranus 

2.3 X 10-11 

2.3 X 10-11 

2.3 X 10-11 

2.4 X 10-11 

2.4 X 10-11 

3.1 X 10-11 

4.5 X 10-11 

2.5 X 10-® 

4.0 X 10-11 

Neptune 

1.7 X 10-11 

1.7 X 10-11 

1.7 X 10-11 

1.8 X 10-11 

1.8 X 10-11 

2.0 X 10-11 

2.5 X 10-11 

4.7 X 10-11 

1.4 X 10-® 


Table 3. The relativistic effect of the oblateness of a body on light at its surface, J2^. Where ^eA s,nd J 2 A ^.re respectively 

c a^A 

the mass, the equatorial radius and the dynamical form factor for body A. 


Sun 

Mercury 

Venus 

Earth 

Mars 

Jupiter 

Saturn 

Uranus 

Neptune 

2.1 X 10-1® 

6.0 X 10-1® 

2.4 X 10-1® 

7.5 X 10-1® 

2.8 X 10-1® 

2.9 X 10-1® 

1.4 X 10-11 

8.4 X 10-12 

5.2 X 10-12 
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Table 4. The order of va between bodies in the Solar System. This table shows that the bodies at columns affect the bodies at rows when 

c^rAB 

a light ray arrives their limbs. 


c-^tab 

Sun 

Mercury 

Venus 

Earth 

Mars 

Jupiter 

Saturn 

Uranus 

Neptune 

Sun 

1.1 X 

10-13 

1.3 X 

10-15 

6.8 X 

10-16 

4.9 X 

10-16 

3.2 X 

10-16 

9.5 X 

10-17 

5.2 X 

10-17 

2.6 X 

10-17 

1.6 X 

10-17 

Mercury 

6.8 X 

10-19 

1.6 X 

10-14 

7.8 X 

10-19 

4.3 X 

10-19 

2.3 X 

10-19 

5.4 X 

10-20 

2.9 X 

10-20 

1.4 X 

10-20 

8.8 X 

10-91 

Venus 

3.9 X 

10-18 

8.4 X 

10-18 

7.0 X 

10-14 

1.0 X 

10-17 

3.5 X 

10-18 

6.3 X 

10-19 

3.2 X 

10-19 

1.5 X 

10-19 

9.6 X 

10-20 

Earth 

2.9 X 

10-18 

4.8 X 

10-18 

1.1 X 

10-17 

6.9 X 

10-14 

5.6 X 

10-18 

7.0 X 

10-19 

3.5 X 

10-19 

1.6 X 

10-19 

1.0 X 

10-19 

Mars 

1.7 X 

10-19 

2.3 X 

10-19 

3.2 X 

10-19 

4.9 X 

10-19 

1.1 X 

10-14 

7.0 X 

10-20 

3.2 X 

10-20 

1.5 X 

10-20 

9.0 X 

10-21 

Jupiter 

7.9 X 

10-1^ 

8.5 X 

10-17 

9.2 X 

10-17 

9.8 X 

10-17 

1.1 X 

10-16 

8.6 X 

10-13 

9.5 X 

10-17 

2.9 X 

10-17 

1.7 X 

10-17 

Saturn 

9.5 X 

10-18 

9.9 X 

10-18 

1.0 X 

10-17 

1.1 X 

10-17 

1.1 X 

10-17 

2.1 X 

10-17 

2.3 X 

10-13 

9.4 X 

10-18 

4.4 X 

10-18 

Uranus 

5.1 X 

10-19 

5.2 X 

10-19 

5.3 X 

10-19 

5.4 X 

10-19 

5.5 X 

10-19 

7.0 X 

10-19 

1.0 X 

10-18 

5.7 X 

10-14 

9.0 X 

10-19 

Neptune 

3.1 X 

10-19 

3.1 X 

10-19 

3.1 X 

10-19 

3.2 X 

10-19 

3.2 X 

10-19 

3.7 X 

10-19 

4.5 X 

10-19 

8.5 X 

10-19 

2.8 X 

10-14 


Table 5. <; and 77 at the spatial isotropic and 

anisotropic parts of c~^ for gij in harmonic gauge 
in different theories. and 7 are PPN parameters 
and ci 4 = Cl + C4, where ci and C4 are constant 
parameters in M theory.^J 


Parameter 

GR 

ST 

M 


1 

272 _ 1 ^ + 2/3 - 1 

1 + ^*^14 
1 — J_C14 

V 

Ref. 

1 

1(1 + 7) 

im 


Table 6. The estimation of the terms in Eq. gsj 


Terml 

Term2 

Term3 Term4 Term5 

Term6 

Term7 

Term8 

Term9 

Term 10 

Term 11 

Term 12 

3.14mas 

18.02/.ias 

1".75 0".47 83.06/ras 

492.76^as 

4.13/ias 

1.53/ias 

3.05^as 

7.44/ras 

11.32//as 

1.20//as 

















